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Abstract 

Electron transport characteristics are investigated through some molecular chains attached to two non- 
superconducting electrodes by the use of Green's function method. Here we do parametric calculations 
based on the tight-binding formulation to characterize the electron transport through such bridge systems. 
The transport properties are significantly influenced by (a) the length of the molecular chain and (b) the 
molecule-to-electrodes coupling strength and here we focus are results in these aspects. In this context 
we also discuss the steady state current fluctuations, the so-called shot noise, which is a consequence of 
the quantization of charge and is not directly available through conductance measurements. 
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1 Introduction 

Recent progress in nanoscience and technology has 
allowed one to study the electron transport through 
molecules in a very controllable way. The prospect 
of individual devices based on single molecules 
is approaching realization since electronic compo- 
nents are getting smaller and smaller. Follow- 
ing experimental developments, theory can play a 
major role in understanding the new mechanisms 
of conductance. The single-molecule electronics 
plays a key role in the design of future nanoelec- 
tronic circuits, but, the goal of developing a re- 
liable molecular-electronics technology is still over 
the horizon and many key problems, such as device 
stability, reproducibility and the control of single- 
molecule transport need to be solved. Electronic 
transport through molecules was first studied the- 
oretically in 1974 pQ. Since then numerous exper- 
iments [21 [3l SJ [6] have carried out on the elec- 
tron transport through molecules placed between 
two metallic electrodes with few nanometer separa- 
tion. It is very essential to control electron conduc- 
tion through such quantum devices and the present 
understanding about it is quite limited. For exam- 
ple, it is not very clear how the molecular transport 
is affected by the structure of the molecule itself 
or by the nature of its coupling with the electrodes. 
The electron conduction through such two-terminal 
devices can be controlled by some bias voltage. The 
current passing across the junction then becomes a 
strongly non-linear function of the applied voltage, 
and its understanding is a highly challenging prob- 
lem. The knowledge of current fluctuations (quan- 
tum origin) provides a key idea for fabrication of 
efficient molecular devices. Blanter et al. [7] have 
described elaborately how the lowest possible noise 
power of the current fluctuations can be determined 
in a two-terminal conductor. The steady state cur- 
rent fluctuations, the so-called shot noise, is a con- 
sequence of the quantization of charge and it can 
be used to obtain information on a system which is 
not available through conductance measurements. 
The noise power of the current fluctuations provides 
an additional important information about the elec- 
tron correlation by calculating the Fano factor (F) 
which directly informs us whether the magnitude of 
the shot noise achieves the Poisson limit (F — 1) or 
the sub-Poisson (F < 1) limit. 

In the present article we provide a simple analyt- 
ical formulation of the transport problem through 
some molecular chains of different lengths by using 
the tight-binding formulation. There exist several 



ab initio methods for the calculation of the conduc- 
tance [El El HOl [111 (HI E] as well as model calcu- 
lations [Tl[I3[Ti[Il[TJ[TJ[2Dl[2Tl[21[13]. The 
model calculations are motivated by the fact that 
the ab initio methods are computationally too ex- 
pensive, while, the model calculations by using the 
tight-binding formulation are computationally very 
cheap and also provide a worth insight to the prob- 
lem. In our present study attention is drawn on 
the qualitative behavior of the physical quantities 
rather than the quantitative ones. Not only that, 
the ab initio theories do not give any new qualita- 
tive behavior for this particular study in which we 
concentrate ourselves. 

The article is specifically organized as follow. Fol- 
lowing the introduction (Section 1), in Section 2 we 
describe very briefly the technique for the calcu- 
lation of the transmission probability (T), current 
(I) and the noise power its current fluctuations (S) 
through a molecular chain sandwiched between two 
metallic electrodes by means of the Green's func- 
tion technique. In Section 3, we study the behavior 
of the conductance as a function of the injecting 
electron energy, the current and the noise power 
of its fluctuations as a function of the applied bias 
voltage for some molecular chains in the aspects of 
(a) the length of the chain and (b) the molecule-to 
electrodes coupling strength. Finally, we draw our 
conclusions in Section 4. 

2 The model and a brief de- 
scription of the theoretical 
formulation 

In this section we describe very briefly about the 
methodology for the calculation of the transmission 
probability (T), conductance (g), current (J) and 
the noise power of current fluctuations (S) for a 
molecular chain sandwiched between two metallic 
electrodes (schematically represented as in Fig. [1]) 
by the use of Green's function technique. 

At low temperature and bias voltage the conduc- 
tance g of the molecular chain is given by the Lan- 
dauer conductance formula [2"4"j . 



where the transmission probability T is expressed 
as [24], 

T = Trace [T s G r M T D G a M ] (2) 

In this expression G r M and G a M are the retarded and 
advanced Green's functions of the molecular chain 
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and Tg and Tjj describe its coupling with the source 
and the drain, respectively. The Green's function 
of the molecular chain is written in this form, 
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where E is the energy of the injecting electron and 
Hm is the Hamiltonian of the chain. This Hamil- 
tonian can be expressed in the tight-binding form 
within the non-interacting picture like, 



H M = H £ic\ci + ^ t 
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where e^'s are the site energies and t is the nearest- 
neighbor hopping strength. In Eq.([3]), Ss and Ed 
correspond to the self-energies due to the coupling 
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Figure 1: (Color online). Schematic view of an ar- 
ray of benzene molecules attached to two metallic 
electrodes (made of gold in the actual experiment), 
the source and the drain, by thiol (sulfur-hydrogen 
i.e., S-H bond) groups (using the chemisorption 
technique) . 

of the molecular chain with the two electrodes. Now 
all the information about the electrode-to-molccule 
coupling are included into these two self-energies 
through the use of Newns- Anderson chemisorption 
theory and the detailed descriptions of this theory 
can be found from the references [HI HI] . 

The current passing across the chain can be as- 
sumed as a single electron scattering process be- 
tween the two reservoirs of charge carriers. The 
current-voltage relationship can be obtained from 
the expression [24 , 



e r 

I(V) = \ \ 
nh J_, 



Us-f D )T{E)dE (5) 



where the Fermi distribution function fglD) = 
/ (E — Hs(D)) with the electrochemical potentials 
^S(D) — Ep ± eV/2. For the sake of simplicity, 
here we assume that the entire voltage is dropped 
across the molecule-electrode interfaces and this as- 
sumption doesn't significantly affect the qualita- 
tive aspects of the current-voltage characteristics. 
This assumption is based on the fact that the elec- 
tric field inside the molecular chain especially for 
short chains, seems to have a minimal effect on the 
conductance-voltage characteristics. On the other 



hand for quite longer chains and higher bias voltage, 
the electric field inside the chain may play a more 
significant role depending on the internal structure 
of the chain [25], yet the effect is too small. 

The noise power of the current fluctuations is cal- 
culated from the following expression [7], 

S = % f°° [T(E){f s (l-f s ) + fn(l-f D )} 
+T(E) {1 - T(E)} (f s - f D f] dE (6) 

where the first two terms of this equation corre- 
spond to the equilibrium noise contribution and the 
last term gives the non-equilibrium or shot noise 
contribution to the power spectrum. By calculating 
the noise power we can determine the Fano factor 
F, which is essential to predict whether the shot 
noise lies in the Poisson or the sub-Poisson limit, 
through the relation [7], 



F 



S 

2^7 



(7) 



The shot noise achieves the Poisson limit for F = 1 
and in this case no electron correlation exists be- 
tween the charge carriers. On the other hand for 
F < 1, the shot noise reaches the sub-Poisson limit 
and it provides the information about the electron 
correlation among the charge carriers. 

Throughout this article we study our results only 
at absolute zero temperature, but the qualitative 
features of all the results are also invariant up to 
some finite temperature (~ 300 K). For simplicity, 
here we take the unit c = e = /i = lin our all 
calculations. 

3 Results and their interpre- 
tation 

Here we will study all the essential features of the 
electron transport characteristics for some molecu- 
lar chains in the two distinct regimes. One is the 
so-called weak-coupling regime, which is defined as 
T s(D) « t and the other one is the so-called strong- 
coupling regime, denoted as ~ t. The symbols 
ts and T£> correspond to the hopping strengths of 
the molecular chain to the source and drain, respec- 
tively. In our present calculations, the parameters 
for these two regimes are chosen as ts = td = 0.5, 
t = 3 (weak-coupling) and t$ — tjj = 2.5, t = 3 
(strong-coupling) . 

In Fig. [2] we plot the variation of the conductance 
g as a function of the injecting electron energy E 
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for the molecular chains in the weak-coupling limit. 
Figures E^a), (b), (c) and (d) correspond to the re- 
sults for the molecular chains with two, three, four 
and five benzene rings, respectively. The sharp res- 
onant peaks are observed in the conductance spec- 
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Figure 2: (Color online). Conductance g as a func- 
tion of the injecting electron energy E for the molec- 
ular bridges in the limit of weak molecular coupling, 
(a), (b), (c) and (d) correspond to the bridges with 
two, three, four and five benzene rings, respectively. 

tra for some particular energy values, while for all 
other energies the peaks almost vanish. At these 
resonances the conductance g approaches to 2 in- 
dicating that the transmission probability T goes 
to unity (since we get the relation g = IT from 
the Landauer conductance formula, see Eq. (JTJ) with 
e = h = 1 in our present formulation). These res- 
onant peaks are associated with the molecular en- 
ergy levels and their positions also change with the 
length of the molecular chain. Accordingly, we get 
more and more resonant peaks as we increase the 
length of the molecular chain. Thus it can be pre- 
dicted that the conductance spectrum manifests it- 
self the electronic structure of the molecular chain. 

In the limit of strong molecular coupling, the res- 
onant peaks in the conductance spectra get substan- 
tial widths as shown in Fig. [3l where Figs.^a), (b), 
(c) and (d) correspond to the same molecular chains 
as in Fig. O Such enhancement of the resonant 



widths is due to the broadening of the molecular en- 
ergy levels caused by the coupling of the molecules 
to the side attached electrodes in this strong cou- 
pling limit, where the contribution comes from the 
imaginary parts of the two self-energies Eg and S^). 
From this figure (Fig. [3|) two important features are 
observed those are: (I) the width and (II) the height 
of the resonant peak across the energy E = grad- 
ually decrease with the increase of the length of 
the molecular chain. Such kind of features are not 
clearly visible in the limit of weak molecular cou- 
pling (see Fig. [2} since the widths of the resonant 
peaks are much narrow for this coupling case. Now 
try to explain these results. With the increase of 
the length of the molecular chain, the two energy 
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Figure 3: (Color online). Conductance g as a func- 
tion of the injecting electron energy E for the molec- 
ular bridges in the limit of strong molecular cou- 
pling, (a), (b), (c) and (d) correspond to the bridges 
with two, three, four and five benzene rings, respec- 
tively. 

levels across the zero energy value (associated with 
the two resonant peaks around the energy E = 0, 
those are observed from Fig. [2]) become closer and 
closer. This is due to the existence of the molecular 
energy levels of the respective molecules at these en- 
ergies and they become broadened with the increase 
of the molecular coupling. Since these two energy 
levels are very closely spaced they overlap with each 
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other in some extent and produce a wide resonant 
peak across the energy E = 0, as shown in Fig. [3] 
The width of this resonant peak decreases gradu- 
ally with the length of the chain which is due to the 
closeness of the two energy levels around the energy 
E = 0. Now the decrease of the height of the reso- 
nant peak across the energy E = is due to the in- 
terference effect of the waves traversing through the 
different arms of the molecular rings. For the chains 
with smaller number of rings, the number of inter- 
fering paths are small, while for the longer chains 
we get several paths. As a result of this interfer- 
ence effect, the probability amplitude of getting the 
electron across the chain becomes strengthened or 
weakened (from the standard interpretation of the 
quantum mechanical theory). The anti-resonances 
in the transmission (conductance) spectra are due 
to the cancellation of the transmittances along the 
two paths. Therefore, the probability amplitude T 
becomes less than one and it gradually decreases 
with the increase of the length of the chain. Such 
anti-resonant states are specific to the interferomet- 
ric nature of the scattering states and do not occur 
in the usual one-dimensional scattering problems in- 
volving potential barriers. From these results it can 
be emphasized that the electron transport charac- 
teristics significantly depend on the length of the 
molecular chain which provides a key idea for fab- 
rication of molecular devices. 

Now we describe the behavior of the current I 
and the noise power of its fluctuations S as a func- 
tion of the applied bias voltage V for these molec- 
ular bridges, where both the current and the noise 
power are evaluated through the integration proce- 
dure of the transmission function T (see Eq.© and 
Eq.©). In Fig. [4) we plot the current and the noise 
power of its fluctuation for the molecular chains in 
the limit of weak-coupling, where Figs. HJa), (b), 
(c) and (d) correspond to the results for the chains 
with two, three, four and five benzene rings, re- 
spectively. The red and blue curves correspond to 
the current and the noise power, respectively. It is 
observed that the current shows staircase-like struc- 
ture with sharp steps as a function of the applied 
bias voltage. This is due to the existence of the 
sharp resonant peaks in the conductance spectra in 
the weak-coupling limit (Fig. [2|) . As the bias volt- 
age increases, the electrochemical potentials on the 
electrodes are gradually shifted and eventually cross 
one of the molecular energy levels. Accordingly, a 
current channel is opened and a jump in the I-V 
curve appears. Another key observation is that for 
the bridge with two benzene rings we get the non- 



zero value of the current after some critical value of 
the applied bias voltage (see the red curve of Fig. 2J 
i.e., the conduction of electron through the bridge is 
started (on-state condition of the molecular bridge) 
beyond some critical value of the bias voltage. This 
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Figure 4: (Color online). Current / (red curve) 
and the noise power of its fluctuations S (blue 
curve) as a function of the applied bias voltage V 
for the molecular bridges in the weak molecule-to- 
electrodes coupling limit, (a), (b), (c) and (d) cor- 
respond to the bridges with two, three, four and five 
benzene rings, respectively. 

threshold voltage of the electron conduction gradu- 
ally decreases with the increase of the length of the 
molecular chain. Thus in the limit of weak molec- 
ular coupling one can tune the threshold bias volt- 
age for electron conduction quite significantly by 
changing the length of the molecular chain. It is 
also noted that the current amplitudes for all such 
bridges are almost invariant since the widths of the 
resonant peaks are nearly same for all these bridges 
(see Fig. [5]) in this coupling limit. 

In the study of the noise power of the current fluc- 
tuations for the weak-coupling limit (blue curves of 
Fig. [4|) we see that for all these molecular wires the 
shot noise goes from the Poisson limit (F = 1) to 
the sub-Poisson limit (F < 1) as long as we cross 
the first step in the current-voltage (I-V) charac- 
teristics. Therefore, we can predict that the elec- 
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trons are correlated after the tunneling process has 
completed. Here the electrons are correlated only 
in the sense that one electron feels the existence of 
the other according to the Pauli exclusion principle 
(since we have neglected all other electron-electron 
interactions in our present formalism) . Another sig- 
nificant observation is that the amplitudes of the 
noise power in the sub-Poisson regime for all such 
wires are almost invariant. These results empha- 
size that in the weak-coupling limit the noise power 
of the current fluctuations within the sub-Poisson 
limit remains in the same level independent of the 
length of the molecular chain. 

The behavior of the current and the noise power 
is also very interesting in the strong-coupling regime 
and in the forthcoming part we will describe about 
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Figure 5: (Color online). Current I (red curve) 
and the noise power of its fluctuations S (blue 
curve) as a function of the applied bias voltage V 
for the molecular bridges in the strong molecule- 
to-electrodes coupling limit, (a), (b), (c) and (d) 
correspond to the bridges with two, three, four and 
five benzene rings, respectively. 

that. In Fig. [5] we display the results of the cur- 
rent and the noise power of its fluctuations in 
the strong-coupling limit for the molecular chains, 
where Figs.[5ja), (b), (c) and (d) correspond to the 
same systems as given in Fig. |4l The red and blue 
curves represent the same meaning as in Fig. [4j The 



shape and height of the current steps depend on the 
width of the resonances in the conductance spec- 
tra. With the increase of molecule-to-electrodes 
coupling strength, the step-like behavior disappears 
and the current becomes a continuous function of 
the applied bias voltage and its magnitude also in- 
creases (see the red curves of Fig. [5J compared to 
the weak-coupling case. Thus one can achieve the 
larger current amplitude by tuning the molecular 
coupling strength, without changing any geometry 
of the molecular bridge. It is also noted that for the 
strong molecular coupling the current amplitude de- 
creases with the increase of the number of benzene 
rings in the chain, which is clearly visible from the 
areas under the curves of the conductance spec- 
tra for the different bridges in this coupling limit 

(Fig. ED. 

Finally, in the study of the noise power of the 
current fluctuations we see that the noise power 
decreases (see the blue curves of Fig. [5]) with the 
number of benzene rings in the chain for this strong 
molecular coupling. The significant observation is 
that for the chains with few number of benzene rings 
(i.e., the chains with two, three and four benzene 
rings), the current fluctuations (see Figs. [SJa), (b) 
and (c)) are in the sub-Poisson limit (F < 1) and 
accordingly, there is no such possibility of transition 
from the Poisson (F — 1) to the sub-Poisson limit, 
like as in the weak-coupling case. Therefore, the 
electron correlations are strongly important for such 
chains (chains with few number of benzene rings) in 
this coupling limit. But for the longer chain i.e., the 
chain with five benzene rings the transition from the 
Poisson limit to the sub-Poisson limit takes place 
beyond the first step in the current-voltage charac- 
teristic (see Fig.[5Jd)). A similar kind of transition 
is also observed for the chains with more than five 
benzene rings in this strong-coupling limit which is 
not shown here in the figure. Thus we can empha- 
size that the noise power of the current fluctuations 
are less significant for the longer chains than that 
of the shorter ones in the limit of strong molecular 
coupling. 

4 Concluding remarks 

In conclusion of this article, we have introduced 
a parametric approach based on the tight-binding 
model to investigate the electron transport charac- 
teristics through some molecular chains attached to 
two metallic electrodes. All the results studied here 
provide several key ideas for fabrication of efficient 
molecular devices. 
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For the weak-coupling limit, the conductance 
shows very sharp resonant peaks (Fig. [2]) associ- 
ated with the molecular energy levels. The widths 
of these resonant peaks get broadened (Fig. [3J in 
the limit of strong molecular coupling, where the 
contribution comes from the imaginary parts of the 
two self-energies Es and Ed [24] , 

In the calculation of the current we have seen 
that the current shows staircase-like behavior (red 
curves in Fig. [4]) with sharp steps as a function 
of the applied bias voltage in the weak-coupling 
limit, while it gets continuous variation (red curves 
in Fig. [5]) with the increase of the molecular cou- 
pling strength. For all such bridges the current am- 
plitudes are comparable to each other in the weak 
molecular coupling, but the current amplitude de- 
cays substantially with the increase of the length of 
the chain in the limit of strong-coupling. 

Finally, in the determination of the noise power 
of the current fluctuations, we have noted that in 
the weak-coupling regime the current fluctuations 
within the sub-Poisson limit are in the same level 
(blue curves in Fig. |4|) independent of the length 
of the molecular chain. On the other hand, for the 
strong-coupling regime the noise power of the cur- 
rent fluctuations are less important for the longer 
chains than that of the shorter ones (blue curves in 
Fig. El). 

Throughout our study, we have ignored the ef- 
fect of electron-electron interaction and the influ- 
ence of all the inelastic processes. More studies 
are expected to take into account the Schottky ef- 
fect which comes from the charge transfer across 
the metal-molecule interfaces, the static Stark ef- 
fect resulting from the modification of the electronic 
structure of the molecular bridge in presence of the 
applied bias voltage (essential especially for higher 
voltages) . 

References 

[1] A. Aviram and M. Ratner, Chem. Phys. Lett. 
29, 277 (1974). 

[2] R. M. Metzger, B. Chen, U. Hoepfner, M. V. 
Lakshmikantham, D. Vuillaume, T. Kawai, X. 
Wu, H. Tachibana, T. V. Hughes, H. Saku- 
rai, J. W. Baldwin, C. Hosch, M. P. Cava, L. 
Brchmer, and G. J. Ashwcll, J. Am. Chem. 
Soc. 119, 10455 (1997). 

[3] C. M. Fischer, M. Burghard, S. Roth, and K. 
V. Klitzing, Appl. Phys. Lett. 66, 3331 (1995). 



[4] J. Chen, M. A. Reed, A. M. Rawlett, and J. 
M. Tour, Science 286, 1550 (1999). 

[5] M. A. Reed, C. Zhou, C. J. Muller, T. P. Bur- 
gin, and J. M. Tour, Science 278, 252 (1997). 

[6] T. Dadosh, Y. Gordin, R. Krahne, I. Khivrich, 
D. Mahalu, V. Frydman, J. Sperling, A. Ya- 
coby, and I. Bar-Joseph, Nature 436, 677 
(2005). 

[7] Y. M. Blanter and M. Biittiker, Phys. Rep. 
336, 1 (2000). 

[8] S. N. Yaliraki, A. E. Roitberg, C. Gonzalez, 
V. Mujica, and M. A. Ratner, J. Chem. Phys. 
Ill, 6997 (1999). 

[9] M. Di Ventra, S. T. Pantelides, and N. D. Lang, 
Phys. Rev. Lett. 84, 979 (2000). 

[10] Y. Xue, S. Datta, and M. A. Ratner, J. Chem. 
Phys. 115, 4292 (2001). 

[11] J. Taylor, H. Gou, and J. Wang, Phys. Rev. B 
63, 245407 (2001). 

[12] P. A. Derosa and J. M. Seminario, J. Phys. 
Chem. B 105, 471 (2001). 

[13] P. S. Damle, A. W. Ghosh, and S. Datta, Phys. 
Rev. B 64, R201403 (2001). 

[14] V. Mujica, M. Kemp, and M. A. Ratner, J. 
Chem. Phys. 101, 6849 (1994). 

[15] V. Mujica, M. Kemp, A. E. Roitberg, and M. 
A. Ratner, J. Chem. Phys. 104, 7296 (1996). 

[16] M. P. Samanta, W. Tian, S. Datta, J. I. Hen- 
derson, and C. P. Kubiak, Phys. Rev. B 53, 
R7626 (1996). 

[17] M. Hjort and S. Staftrom, Phys. Rev. B 62, 
5245 (2000). 

[18] K. Walczak, Phys. Stat. Sol. (b) 241, 2555 
(2004). 

[19] K. Walczak, Cent. Eur. J. Chem. 2, 524 (2004). 

[20] D. Walter, D. Neuhauser, and R. Baer, Chem. 
Phys. 299, 139 (2004). 

[21] R. Baer and D. Neuhauser, Chem. Phys. 281, 
353 (2002). 

[22] R. Baer and D. Neuhauser, J. Am. Chem. Soc. 
124, 4200 (2002). 



7 



[23] K. Tagami, L. Wang, and M. Tsukada, Nano 
Lett. 4, 209 (2004). 

[24] S. Datta, Electronic transport in mesoscopic 
systems, Cambridge University Press, Cam- 
bridge (1997). 

[25] W. Tian, S. Datta, S. Hong, R. Rcifcnberger, 
J. I. Henderson, and C. I. Kubiak, J. Chcm. 
Phys. 109, 2874 (1998). 



8 



